A note on bound entanglement and local realism. 
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Abstract: We show using a numerical approach that gives necessary and sufficient conditions for the existence of 
local realism, that the bound entangled state presented in Bennett et. al. Phys. Rev. Lett. 82, 5385 (1999) admits a 
local and realistic description. We also find the lowest possible amount of some appropriate entangled state that must 
be ad-mixed to the bound entangled state so that the resulting density operator has no local and realistic description 
and as such can be useful in quantum communication and quantum computation. 

r— I ■ I. INTRODUCTION 

o: 

*vj Since the discovery of the protocol of quantum teleportation in 1993 lit] it has been shown that quantum entangle- 

^^ ment is a fundamental resource in quantum communication. Some recent research suggests that the phenomenon of 

^ , quantum entanglement can be also viewed as a resource in quantum computation H. 

' ' The best source of quantum entanglement is, of course, a pure maximally entangled state. However, if one thinks 

^-H _ about quantum communication and quantum computation in a noisy environment than it is necessary to deal with 
f^ ■ mixed states. It has been shown that for a wide class of mixed entangled states by distillation protocol H one can 
, , extract useful entanglement that can be then used as a resource in quantum communication and quantum computation. 
>^ • Nevertheless, there still exists a set of entangled mixed states which cannot be distilled B . This phenomenon has 
T^ I been called bound entanglement. Although, these states are undistillable they are not completely useless in quantum 
lO ■ communication. In |5[| it has been shown that some family of bound entangled states can be useful in conclusive 
quantum teleportation process. 

The discovery of bound entangled states has raised a question whether they can be described in a local and 
realistic way, i.e., if there exists a Bell- type experiment in which predictions of local realism are violated by quantum 
^ ■ mechanical ones. It has been Peres who first conjectured that such states should admit a local and realistic description 
^ Q. However, as there are no known Bell inequalities being a sufficient and necessary condition for the existence of 

(-^ local realism for systems described by the tensor product of two Hilbert spaces each of the dimension greater than 
I , two (we know that lowest dimensional bound entangled state lives in a tensor product of two Hilbert spaces of the 
(~| ' dimension equal to three), it could not be provedg. Therefore, the question has remained open. 

Ci ' In this paper we show that using numerical method of linear optimization we can answer the above question. 

Although, here we investigate the specific bound entangled state the presented method can be successfully applied 
to any bound entangled state (in general, to any mixed state). The strength of the method lies in the fact that it 
gives necessary and sufficient conditions for the existence of local and realistic description of the investigated quantum 
system. 
5^] ' Additionally we calculate how much of some optimally defined entangled state must be ad-mixed to the investigated 

Cu I bound entangled state so that the created mixture does not have a local and realistic description anymore. 

II. CONSTRUCTION OF THE BOUND ENTANGLED STATE 

In ra] it has been shown a simple and elegant way to construct bound entangled states. The idea is based on the 
theorem proved by Horodecki M and on the notion of unextendible product basis (UPB) m. 

In short, the Horodecki theorem states that if a density matrix p of two quantum systems A and B is separable 
then there exists a family of product vectors, say \(f>k)\'>Pk) (first ket refers to system A and the second one to system 
B), spanning the range of p and such that vectors |0fc)|'0fc)* span the range of partially transposed matrix p'^^ , where 
Tb means the transposition with respect to the subsystem B and the star denotes complex conjugation. 
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'^Recently discovered Bell inequalities for higher dimensional quantum systems [y,|lO| have not been proved to be sufficient 
condition for the existence of local realism 



An unextendible product basis is a set of d orthogonal product vectors belonging to A'' x M dimensional Hilbert 
space, where d < N x M, and of the property that any orthogonal vector to them must be entangled. 

Using the definition of UPB and the Horodecki theorem one can easily construct bound entangled states. Let us 
assume that we have a Hilbert space being a tensor product of two N dimensional Hilbert spaces. In this space we 
find some unextendible product basis consisting of d product vectors |0fc)|'0fc), where fc = 1, 2, . . . , d. Denoting by P 
the projector on the subspace spanned by them one can define the following density matrix p 

P = T^Ai-n (1) 

The support of this matrix (space spanned by its eigenvectors) lies in a subspace perpendicular to the subspace 
spanned by the vectors from the unextendible product basis, therefore it is spanned by -/V^ — d entangled orthonormal 
vectors. The Horodecki theorem says that this state must be entangled and as it can be directly checked its partial 
transposition is a positive operator. Hence we have obtained a bound entangled density matrix. 

In this paper we investigate a particular bound entangled state consisting of two subsystems each living in a three 
dimensional Hilbert space. As the UPB wc choose the following set of orthonormal vectors ^ 

\v,) = ±i\0)\0)-\0)\l)) 

|z;i) = i=(|0)|2)-|l)|2)) 

|z;2) = -^(|2)|l)-|2)|2)) 

|z;3) = (^(|l)|0)-|2)|0)) 

|z;4) = ^(|0) + |l> + |2»(|0) + |l> + |2>, (2) 

and follow the recipe described above. The resulting density matrix can be written as follows 

PB = -^{\V5){V5\ + MM + \V7){V7\ + \vs){v8\), (3) 



where 



1^5) = -^(K) \vt)) 



1 

V2' 



\ve) = —{\V + ) - ^3+)) 



\v7) = l{H) + \vt)-\vt)-\vt)) 

\vs) = ^(1^0+) + 1^1+) + \4) + H)) ^|1>|1), (4) 

and where vectors |w^) (fc = 0, 1, ... ,3) are made of the states belonging to (||) by changing the sign — to +. For 
instance, \v+) = -i=(|0)|0) + |0)|1) and so on. 

The above vectors are normalized, orthogonal to each other and orthogonal to vectors forming the UPB basis, i.e., 
vectors \vi) with / = 0, 1, 2, . . . , 4. Therefore, they span the orthogonal complement of the subspace associated with 
the UPB basis. As a whole, the vectors \vk) with fc = 0, 1, . . . , 8 form an orthonormal basis in Hilbert space. 

III. LOCAL HIDDEN VARIABLES 

On the above mixed state ps one can perform a Bell type experiment in which two spatially separated observers 
Alice and Bob measure some trichotomic observables. In this paper we consider the case in which both observers 
are allowed to measure only two noncommuting trichotomic observables, which we denote by Ai,A2 for Alice and by 
Bi,B2 for Bob. 



Any trichotomic observable in a three dimensional Hilbert space can be obtained by a rotation of some orthogonal 
basis by means of a unitary transformation belonging to SU{3) group. SU{3) group is a set of unitary matrices with 
the determinant equal to one and which depend on eight real parameters 4'i,4>2, ■ ■ ■ ,4>i,, which we will denote as a 
vector (j) = ((/)i, (J)2t ■ ■ ,4's,)- Suppose that the resolution of unity at Alice's side consists of three orthogonal projectors 
Pfc = |^)(^| and that at Bob's side of projectors Qi — \l){l\ (fc, Z = 1, 2, 3). Then arbitrary trichotomic observables for 
Alice and Bob read 

i. = alC/(4)Ac/(4)t + a\U{$\)p2U{$\)^ + c^^U{$\)P:iU{$\)^ 

B, = 6iy(4)A^(0l)^ + 6^2^(^i)AF(<?i)t + biV{$^s)P:,V{$'s)\ (5) 

where C/((/)^), y((/)-^) (i,j = 1,2) are members of SU{3) group and where the numbers a^,6|, {k = 1,2,3) are the 
eigenvalues of appropriate observables. Of course each observable Ai and Bj depends on vectors 4>\ and (j)-'^ but to 
shorten notation we will not write it explicitly. 

The probability PQM{ak,bj\4>\,(j>''g) of obtaining the pair of eigenvalues a^,6f while measuring the observables 
Ai , Bj on the density matrix p can be calculated in a standard way as 

PQMiam\$\J%) = TriUi$\)PkUi$\y®Vi$^)QiVi$%yp). (6) 

Obviously, there are 36 such probabilities. 

A local and realistic description of the presented quantum experiment is equivalent to the existence of a joint 
probability distribution, which returns all 36 quantum probabilities as marginals [|ll| . Let us denote this hypothetical 
probability distribution as PLfl(a^, a^, 6^, 6^), where k,l,m,n — 1,2,3. It consists of 81 non-negative numbers 
summing up to one. The marginals are given by the following set of equations 



PLRialbl) = J2J2 PLR{al,albl,bl) 

I n 

PLRialbl) ^Y.T. PLR{al,albl,bl) 

I m 

PLR{albl,) ^Y.11 PLR{alalbl,bl) 

k n 



k rn 



A quantum experiment admits a local and realistic description if and only if the above marginals can be made equal 
to quantum ones ||ll|, i.e., PLR{al,bl) = Pqm«, 6f |(^^, ^■^). 

As there are not known Bell inequalities for two three dimensional quantum systems that are necessary and sufficient 
condition for the existence of local realism one is forced to resort to numerical methods. The problem can be solved 
numerically by means of linear programming method |12|Jl3|| . Linear programming is a method of finding the maximum 
of a linear multi- variable function, which domain is a convex set. It relies on the fact that the maximum is reached 
(if it exists) in one of the vertices of the domain. Moreover, linear programming can be used (by setting the function 
being maximized to be constant) to check if the set of linear equations that form the boundaries of a convex set 
(domain) has a solution. In this case the lack of the solution means that the considered convex set is empty, i.e., it 
has no " interior" . 

Therefore, linear programming can be successfully applied to our case. Indeed, we are interested in finding a local 
realistic probability distribution PLfl(a^,a^,6^, 6^), i.e., 81 non-negative numbers, returning quantum probabilities 
as marginals. The set of equations (0) plus the condition that the probabilities sum up to one defines, for the given 
set of observables, a convex set in an 82 dimensional real vector space. Now, applying linear programming procedure 
we can check if the set is empty. If not, then there exists a local and realistic description of the experiment. 

Thus, for each choice of the observables Ai, A2, Bi, B2 we apply the numerical linear programming procedure. In 
our calculations we have used the state-of-the-art HOPDM 2.30 procedure [Ea. The quadruples of observables (two 
for Alice and two for Bob) have been chosen randomly 10^ times. Local and realistic description has existed for all 
these cases. Therefore, it suggests that the entanglement contained in this state is too weak to violate local realism 
and as such it is a strong evidence supporting Peres' conjecture [pj that bound entangled states admits a local and 
realistic description. 



IV. EXTRACTING ENTANGLEMENT 

Although the number (10®) of Bell experiments simulated numerically is huge it may be that one can still find some 
quadruple of observables for which there is a violation of local realism by the bound entangled state pB- Therefore, 
it is desirable to support the above result using some additional calculations. To this end we propose the following 
procedure. 

We ask what is the least amount F (0 < F < 1) of some entangled states |'0)(V'| that has to be admixed to the 
state pb so that the resulting state p{F) = (1 — F)pB + F\il)){ip\ has not a local realistic description anymore. It is 
clear that in such a case our numerical procedure for large F cannot return a local hidden variable model for any set 
of local settings. However, for lower _F's, as it turns out, it does find such models. Since numerical linear optimization 
is a highly reliable procedure, and we additionally use another procedure (see below) to find the optimal parameters 
of the problem for the existence of local hidden variable models, the results obtained in this way are much more 
definitive than the "lottery" approach presented above. Also they in a way measure the "robustness" of the local 
realistic models for the bound entangled state. 

As the support of the state pB is spanned by the vectors \vi) [i — 5,6,7,8) it is naturally to consider the 
state \ip) as a superposition of these vectors, i.e., \ip) — X]i=5 '^il^i)j where X]i=5 I'^iP = 1- It is convenient to 
parametrize the complex numbers a^ by six angles ^,^, (/>, Xi: X2j X3 i^ the following way: 05 ~ sin f/; sin cos (/i, ag — 
exp(zxi) sin?/;sin0sin(/), 07 = exp(ix2) sinV'cos^, ag = exp(ix3) cos V'- 

To find the optimal state jV') (optimal in the sense defined above) for every choice of angles "^j^j </'jX1: X2: X3 
and the observables Ai,A2,Bi,B2 (we remind that each observable depends on 8 angles) we calculate using linear 
programming procedure HOPDM 2.30 the maximal value of the parameter F, which now depends on 38 angles 
(32 angles defining observables and 6 angles defining the state \ip)), below which there exists a local and realistic 
description of the experiment. This way we obtain the 38 variable function which minimum Fmin can be found by 
the so called amoeba procedure utilising downhil simplex method ||l4| . Although the amoeba is a good minimization 
procedure there is no guarantee that the found minimum is a global one (a problem encountered in any numerical 
minimization). To reduce the risk of finding the local minimum the procedure has been run many times with randomly 
chosen initial conditions. 

Calculations show that the minimal possible F equals F„iin — 0.509651. This occurs for the state \^p) — -y^{\v5) + 
\ve)) (the analytical form of the state [ip) has been obtained using the numerical results and then, for additional 
confirmation, the numerical calculations have been performed again for the guessed state without the minimization 
over the angles characterizing it yielding the same Fmin)- 

It can be directly checked that the state p{F) for any F 7^ is negative under partial transposition. Therefore, we 
have the situation in which the family of mixed states p{F) for < F < Fmin is entangled but has a local realistic 
description. This once again shows that non-separability is not equivalent to lack of classical description. 

It is instructive to calculate the degree of entanglement of the state jV') defined as the |(1 — Frs((Fr^(|?/;) (V'l))^)), 
where for instance TrA denotes the trace with respect to Hilbert space of Alice's subsystem. Trivial algebra gives 
us y|. It can be checked numerically that this is the most entangled state that can be obtained by superposing the 
states \vi) (i = 5,6,7,8). This is an interesting confirmation of something that one could intuitively expect. The 
most efficient way to get to a region in which there is no local realistic description is via an admixture of the most 
entangled pure state that lives in the subspace of the Hilbert space which is used to construct the considered bound 
entangled state. 

V. CONCLUSIONS 

The state ps cannot be distilled, which means that the entanglement is hidden in the state too deep to be recovered 
by local quantum operations. Moreover, the presented here numerical computations strongly support the hypothesis 
that the bound entangled state pB is describable by local hidden variables or in other words that the correlations 
observed in a Bell-type experiment with this state can be simulated classically. The large value of F, the admixutre 
of optimal entangled state required to get violations of local realism, clearly indicates that the considered bound state 
is not even close to a border of the realm of states possessing a local realistic model. It is well inside this realm. 

It has been recently shown that multipartite (more than seven quantum systems) bound entangled states violate 
local realism |L6| . In the view of the results presented here, i.e., no violation of local realism by the bipartite bound 
entangled state pB , one sees that one is far away from a complete understanding of the structure of bound entangle- 
ment. We hope that it will stimulate further research in this direction. 
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